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ABSTRACT 

Context. The Ly a line is an important diagnostic for star formation at high redshift, but interpreting its flux and line profile is difficult 
because of the resonance nature of Ly a. Trends between the escape of Ly a photons and dust and properties of the interstellar medium 
(ISM) have been found, but detailed comparisons between Ly a emission and the properties of the gas in local high-redshift analogs 
are vital for understanding the relation between Ly a emission and galaxy properties. 

Aims. For the first time, we can directly infer the properties of the ionized gas at the same location and similar spatial scales of the 
extended Ly a halo around the local Ly a emitter and Lyman-break galaxy analog ESO 338-IG04. 

Methods. We obtained VLT/MUSE integral field spectra. We used ionization parameter mapping of the [S n]/[0 hi] line ratio and the 
kinematics of Her to study the ionization state and kinematics of the ISM of ESO 338-IG04. 

Results. The velocity map reveals two outflows, one toward the north, the other toward the south of ESO 338. The ionization parameter 
mapping shows that the entire central area of the galaxy is highly ionized by photons leaking from the H n regions around the youngest 
star clusters. Three highly ionized cones have been identified, of which one is associated with an outflow detected in the Her. We 
propose a scenario where the outflows are created by mechanical feedback of the older clusters, while the highly ionized gas is caused 
by the hard ionizing photons emitted by the youngest clusters. A comparison with the Ly a map shows that the (approximately bipolar) 
asymmetries observed in the Lycr emission are consistent with the base of the outflows detected in Her. No clear correlation with the 
ionization cones is found. 

Conclusions. The mechanical and ionization feedback of star clusters significantly changes the state of the ISM by creating ionized 
cones and outflows. The comparison with Lycr suggests that especially the outflows could facilitate the escape of Lycr photons. 
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1. Introduction 


Lycr is one of the most important diagnostic lines in extra- 
galactic astrophysics. The line is thought to be bright in star¬ 
forming galaxies, allowing their detection at high redshift. How¬ 
ever, due to the resonance nature of the Lycr line, relating the line 
strength to star formation is not a straightforward process, and 
Lycr needs to be used only with great care as a physical probe. 
A Lycr photon scatters in neutral H i gas until it escapes or is ab¬ 
sorbed by dust. This means that Lycr is a complicated function 
of not o nly the ionizing power and extincti on, but also of the 
amount (iPardv et al1l2014l) and kinematics dKunth et al.l Il998b 
Woff ord etal.1120131) of H i. 


To understand the relation between Lycr and these properties, 
detailed studies of local analogs of Lycr-emi tting galaxies are 
needed (lOstlin et al.ll20l4 lHaves et aDl2014l) . Observations of 
nearby galaxies have shown that stellar feedback creates highly 
ionized cones in the neutral medium, allowing the ionized radia¬ 
tion to escape from the galaxy (Zast row et al.fl20 1 ll 120 1 3l) . This 
would make the amount of escaping Lyman-continuum (LyC) 
photons a strong function of the orientation under which we ob¬ 
serve the galaxy. 


* Based on observations collected at the European Southern Obser¬ 
vatory at Paranal, Chile (ESO program 60.A-9314). 


In this letter, we compare for the first time the spatially re¬ 
solved Lycr emission, obtained with th e Hubble Telescope (HST) 
(lHaves et al.ll2005l : lOstlin et al.ll2009h . with the ionization state 
and kinematics of the ionized gas in the halo of the bright¬ 
est Lycr emitter in the local Universe and Lyman-break galaxy 
(LBG) analog, the starburst galaxy ESO 338-IG04 (Tololo 1924- 
416, hereafter ESO 338). We present high-quality MUSE ob¬ 
servations allowing a detailed comparison between the spatial 
distribution of the Lycr emission and potential outflows in this 
galaxy. In Sect. 2 we describe the observations and data reduc¬ 
tion. In Sect. 3 we present the results of the ionization mapping 
and kinematics of Her toward the galaxy. We conclude with dis¬ 
cussion and conclusions in Sect. 4. 

2. Observations and data reduction 

ESO 338 has been observed with the integral field spectro grap h 
Multi-Unit Spectroscopic Explorer (MUSE) dBacon et aL]l2010l) 
mounted on the Very Large Telescope (VLT) on Paranal, Chile. 
The observations were performed during the first science verifi¬ 
cation run. The observations are taken in the no-AO mode in the 
extended wavelength setting, providing spectra between 4600 
and 9350 A in a field of view of 1 □' with a 0.2"pixel scale. 
The observations were performed on 2014 June 25, split over 
two observing blocks with four frames of 750 sec, rotated by 
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Fig. 1 . Left: Reconstructed VRHcr color image of ES0338 as seen by MUSE. The Her image is not continuum subtracted. The black asterisks 
mark the location of the two foreground stars. The white box shows the region displayed in the panel to the right. Right: Reconstructed three-color 
images based on narrow-band images taken with HST of Lya (blue), [O hi] 25007 A (green) and Her (red). 


90°, resulting in an integration time of 100 min. The measured 
seeing was 0.9"in the V band. The data were reduced using the 
ESO pipeline version 0.18.5 (Weilbacher et al, in prep.). Flux 
calibration is based on observations of the spectrophotometric 
standard star GD 153. 

We present the first results derived from this rich data set. 
Broadband images in V and R were extracted using the V and 
R transmission curves. Emission line maps of [Sn] 26717 A 
and 26731 A, [Om] A5007 A, and H a 26563 A were extracted. 
The emission line maps were constructed by numerically inte¬ 
grating below the emission line, and the continuum was sub¬ 
tracted by averaging the continuum blue- and redward of the 
line. To enhance the low-surface brightness features, we made 
use o_f_the weighted Voronoi tesselations binning algorithm by 
Diehl & Statler ( 2006) , which is a gene r alizati on of the algo¬ 
rithm developed bv ICappellari & Copinl (120031) . A minimum 
signal-to-noise ratio (S/N) of 20 has been chosen for the [S n] 
26717 A line with a maximum cell area of 100 pixels (4 □"). 
This Voronoi pattern was applied to all emission line maps. A 
Gaussian profile was fitted to the H a line using the MPFIT 
(lMarkwardl2009l) to calculate the velocity map. 


3. Results 

Figure \T\ (left) shows a VRHa three-color image extracted from 
the MUSE observations of ESO 338. The reconstructed V and R 
images reveal the continuum emission of the galaxy and show 
numerous star clusters (lOstlin et ai1ll998l). The H a map show s 
a large ionised halo around the galaxy dBergvall & Ostlin 2002), 
extending as far as 6 kpc from the center of the galaxy (assuming 
a distance of 37.5 Mpc. ilOstlin et al. 19981 . 

The right panel of Fig. |T] shows a composite of HST obser¬ 
vations of H a, O hi, and Ly a of the central part of ESO 338. The 
MUSE observations allow us to study the gas far out in the halo, 


while the HS T observations give detailed informa tion about the 
small scales (lHaves et al.ll2005l: lOstlin et al.ll2009l) . 

3.1. Ionization parameter 

We applied the technique of ionization parameter mapping to 
determine the degree of ionization and the optical depth for Ly a 
photons using the [SII] (26717 A+26731 A)/[OIII] (25007 A) 
line ratio dPellegrini et al.1120121) . The right panel of Fig. [2| shows 
the [S n]/[0 hi] line ratio calculated using the Voronoi tessellated 
data. The observed ratio varies from unity in the western and 
eastern parts of the halo to 0.03 in the center of ESO 338. This 
extreme line ratio shows that there is a high production of hv > 
34 eV photons compared to where [S n]/[0 hi] is L Comparing 
these low values with theoretical prediction by [Pellegrini et al.1 
( 2012 ) shows similarities to the models with an 03V star as 
ionizing star. The lowest values are found coincident with the 
positions of young (<2 Myr) star clusters (lOstlin et al.ll2003l) . 
containing these very early O stars. 

The central ~2 kpc of the galaxy halo show an average value 
of 0.1 of the [Sii]/[Oiii] ratio, indicating that entire center of 
ESO338 is highly ionized. This is probably caused by the Hu 
regions being optically thin for LyC photons, which allows them 
to leak out and ionize of the gas outside the H n regions. The gas 
becomes optically thick for LyC photons only in the western and 
eastern outer regions of the galaxy. This pattern mimics the be¬ 
havior of a gian t H n region, as found for similar galaxies (Haro 
11. ICormier et al.f2012 ). and confirms the existence of an in-out 
ionizing front that modifies the interstellar medium (ISM) states 
of the galaxy. Several cones of highly ionized gas expanding all 
the way into the galaxy halo are visible. The most prominent ion¬ 
ization cones are in the northwest (nr 1) and southeast direction 
(nr 3). An additional ionization cone (nr 2) is found in the eastern 
direction. 
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A. Bik et al.: VLT/MUSE view of the highly ionized outflow cones in the nearby starburst ESO338-IG04 
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Fig. 2. Left: Velocity map of Her as derived from a Gaussian fit of the emission line. As systemic velocity we have chosen 2841 km s _1 . Overplotted 
are the contours of the Lya -continuum subtracted emission (lOstlin et al.ll2009f) w ith con tour levels of 0.3, 0.9 and 2 x 10 -19 erg/s/cm 2 /A. Right: 
[S n]/[0 hi] line ratio tracing the ionization parameter throughout the galaxy ([Pellegrini et al. 11 1201 2). All cells with a flux lower than 1.7 x 10“ 19 
erg/s/cm 2 /A, corresponding to a S/N = 20 for the S ii 46717 A line, have been removed from the plot. The contours show the Her velocity map 
with contours at -15 (dashed), 15, 30, and 45 (solid) km s -1 ). 




Changes for instance in extinction and abundance might also 
result in changes in line ratios. However, the extinction toward 
ESO 338 is very low, with a luminosity-averaged value of close 
to 0 mag (lOstlin et al.1120031) . This is confirmed by our own 
H/3/H a map, where maximum values of E(B - V) = 0.2 mag 
are measured toward some clusters. To reproduce the observed 
change in [S n]/[Om] ratio, E(B - V ) values as high as 3.5 mag 
are needed dCalzetti et alJll 2000 h . To test whether the abundance 
variations might cause the observed pattern, we calculated the 
N 2 and S 3 O 3 indices (IStasinskal 112006 ). Spatial variations were 
found, but not in the same geometry as for [S n]/[0 111 ]. The 
same spatial variations, however, were found in line ratio maps 
of lines originating from different ionization potentials, such as 
H(r/[Om] and [S n]/[S 111 ], suggesting that the changes in the 
[S n]/[0 hi] ratio are dominated by variations in ionization. 


3.2. Kinematics 


To relate the detected ionization cones to galactic-scale outflows, 
we compared the map of ionization parameters with the H a ve¬ 
locity map. Previous studie s on the Hot velocity fiel d of ESO 338 
based on Fabr v-Perot data (lOstlin et al.ll999[l200H) and long-slit 
spectroscopy dCumming et al.l 120081) concluded that either the 
galaxy is not in equilibrium, or the H a velocities do not trace 
the gravitational potential and are dominated by feedback. 

With our high-quality MUSE data we clearly detected sev¬ 
eral high-velocity outflows (Fig. [2]) and confirm that most of 
their observed velocity profile is indeed dominated by the out¬ 
flow emission. As a result of the irregular velocity profile, the 
systemic velocity is only poorly defined. We measured the ve¬ 
locity of the H a line of the intensity-weighted average over 
the part of the galaxy where continuum is visible. This re¬ 
sulted in a observed velocity of 2841 ± 1 km s -1 , consistent 
with HIPASS 21 cm measurements ( Meyer et al . 1200 4). Mea¬ 
surements of absorption lines of Nai (4 5889.9A, 5895.9A) and 
Mgi (4 5167.3A, 5172.7A, 5183.6 A) in the intensity-weighted 


spectrum show similar values, albeit at lower accuracy because 
of their low equivalent width. This is consistent with previous 
studies where no velocity difference was found between stel¬ 
lar absorption lines and neb ular emission lines ( Cumming et al.l 
12008k ISandberg et al.l l2013h . After correcting for the heliocen¬ 
tric velocity, we chose the systemic velocity to be 2841 km s -1 
as measured on the H a line of the integrated spectrum (Fig.0. 

A complex of outflows coincides with the northwestern ion¬ 
ization cone, suggesting that this cone is created by the expand¬ 
ing, highly ionized gas that is expelled from the galaxy by stellar 
feedback. A similar outflow is found to the south of the galaxy, 
but this outflow is not aligned with an ionization cone. In ad¬ 
dition, two expanding arcs are detected in the velocity map, 
one redshifted, roughly coinciding with ionization cone 3, and 
a blueshifted arc west of the same ionization cone. The out¬ 
flows show maximum projected velocities of 50 km s^ 1 , which 
is slo wer than outflo ws in M82 (IShopbell & Bland-H aw thorn 
1998| ) and NGC 1569 dWestmoquette et al.ll2008l) . for instance, 
which were derived from data with much higher spectral reso¬ 
lution. We found no evidence for multiple components in H a, 
probably because of the insufficient spectra resolution. 

The MUSE data are much deeper than the previous spectro¬ 
scopic H a observation s of ESO 338 of lOstlin et al.l (Il999h and 
ICumming et al.l d2008). This allowed us to trace the halo much 
farther out than before (6 kpc, V- 11 kpc). In addition to the out¬ 
flows, we find evidence for a gradient in the velocity of the halo 
gas. The observed velocities increase from west to east between 
-50 km s -1 and +10 km s -1 , similar to the gradient observed in 
the H 1 emission of ESO 338 (ICannon et al.ll2004l) . 

4. Discussion and conclusions 

We presented MUSE integral field observations of the local Lya- 
emitting galaxy ESO 338. We identified several highly ionized 
channels and found evidence of galaxy-scale outflows driven by 
stellar feedback in the central part of the galaxy. ESO 338 is per- 
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fectly suited for studying the relation between the Ly a escape 
and ionization channels and outflows because it is the LBG ana¬ 
log closest to us. 

4.1. Feedback 

A comparison between the ionization cones and the outflows 
shows that there is no one-to-one match. The northern outflow 
aligns with ionization cone nr 1. The southern part of the galaxy 
halo shows a different geometry. The outflow and ionization 
cone nr 3 are not aligned and are probably two unrelated features. 
We did not detect velocity changes for ionization cone nr 2. 

The two arc features detected in the velocity map might be 
part of the same expanding structure. Figure [2] shows that the 
blueshifted arc marks the edge of the highly ionized region, 
while the red arc is located at the back side of the galaxy ex¬ 
panding away from us. With this hypothesis, they could trace 
the shock front created by the expansion of the hot, overpres¬ 
sured gas from the center of the galaxy into the more neutral 
ISM. 

Massive stars are sources of feedback, therefore they are the 
main agents of the modification of the ISM on both local and 
galactic scales. The large majority of massive stars are in clusters 
because ESO 338 has a very high cluster formation efficiency 
150%. lAdamo et alJil2Ql 1 ). Hard ionizing photons are only pro¬ 
duced by the most massive stars, which in turn means that they 
are produced only in clusters with ages younger than ~ 4 Myr, 
while the mechanical energy released by stellar winds and super- 
novae remains roughly co nstant for ~ 30 Myr dLeitherer et al.l 
1 19991) . lOstlin et all (120031) derived the star formation history 
based on the detected clusters and showed that the star forma¬ 
tion increased 20-30 Myr ago and again 6-10 Myr ago, and it 
increased a little during the last 3 Myr. The oldest star forma¬ 
tion burst has injected the most mechanical energy into the ISM 
and therefore might be responsible for several of the observed 
outflows. The LyC photons responsible for the ionized cones 
are created by the most recent increase in star formation. These 
channels might be outflow cones created by the older generation 
of clusters and might be re-ionized by the younger generation. 

That not all ionization cones are associated with outflows 
might also be caused by projection effects. The ionization cones 
are detected at any inclination because what we measure is the 
average ionization along the line of sight. Only if the cone origi¬ 
nates in the center and is directed straight toward us will the cone 
not be detected because the center is already highly ionized. The 
detection bias introduced by velocity measurements is different, 
such that if the outflows are in the plane of the sky, they will 
not be detected in the velocity map. A more detailed comparison 
with higher spatial resolution between the outflows and clusters 
is needed to test these hypotheses. 

4.2. Lya escape 

The Lya emission map of ESO 338 (Fig. |T]) shows strong emis¬ 
sion and absorption toward the young star clusters in the cen¬ 
ter of the galaxy dHaves et al .1120051) . The Lya photons are pro¬ 
duced in the H n regions surrounding the star clusters. The low 
[S n]/[0 hi] ratio of the central area of the entire galaxy shows 
that these H n regions are optically thin and leak LyC photons. 

As noted bv lHaves et al. ( 2005 ). the clusters in ES0338 are 
surrounded by a diffuse halo of Lya caused by the scattering 
of Lya photons in the neutral ISM. In addition to this scatter¬ 
ing halo, two flux enhancements in the Lya map are observed as 


well: One at the northern side of the galaxy and one toward the 
southwe st, and a relation with possible outfl ows has been pro- 
posed (lHaves et al.1120051: lOstlin et al.ll2009h . With the MUSE 
data we can relate these asymmetries to the outflows detected in 
the H a velocity map. The inner parts of the two detected out¬ 
flows in the velocity map line up with the Lya asymmetries. The 
expanding gas in the outflow causes the resonance wavelength 
to shift away from the Lya wavelength, which allows the Lya 
-photons to escape instead of being scattered. 

The observed outflows as well as the geometry of the Lya 
emission shows that the determination of the Lya escape fraction 
will depend strongly on the ISM geometry of the galaxy. If we 
observe the galaxy through an outflow cone, the observed Lya 
flux may be much higher than for the way we observe ESO 338. 
This allows testing theoretical models that predict that the varia¬ 
tions o f Lya are a function of the viewing angle of the ionization 
cones (lLaursen et al.ll2009l:lBehrens et al.ll2014l) . 

We did not observe a correlation with the ionization cones, 
apart from where the ionization cone aligns with the northern 
outflow. This might suggest that outflows are more effective in 
enhancing the Lya escape fraction than ionization cones, which 
are more likely to facilitate the LyC continuum escape (e.g., 
IZastrow et al.ll2013l) . Interestingly, ESO 338 is also a good can¬ 
didate galaxy in which to study the emission of LyC radiation. 
iLeitet et al . 3201 31) estimated a LyC escape fraction of -16% 
based upon the residual intensity in the usually saturated Cn 
A = 1036 A absorption line. 
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